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Abstract 



^..^ , An amplitude equation for an unstable mode in a collisionless plasma is 

r— ^ \ derived from the dynamics on the two-dimensional unstable manifold of the 

o, 

■^ \ equilibrium. The mode amplitude p{t) decouples from the phase due to the 



0\ 
O 

C/3 



spatial homogeneity of the equilibrium, and the resulting one-dimensional 

dynamics is analyzed using an expansion in p. As the linear growth rate 7 
-(— > , 

o ' vanishes, the expansion coefficients diverge; a rescaling p{t) = 7^ ^(7^) of the 



mode amplitude absorbs these singularities and reveals that the mode electric 
J^ I field exhibits trapping scaling |ii^i| ~ 7-^ as 7 — > 0. The dynamics for r{T) 

depends only on the phase e*^ where dck/dz = |efc|e~*^'^ is the derivative of 
the dielectric as 7 — > 0. 
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The collisionless evolution of an unstable electrostatic mode is a fundamental topic in 
the theory of strong wave-particle interactions. The linear instability arises from a resonant 
interaction between an initial electrostatic fluctuation and particles at the phase velocity of 
the linear mode; the nonlinear evolution is marked by the trapping of the resonant particles 
in the wave potential and decaying oscillations in the wave amplitude due to the bouncing 
and phase mixing of the trapped particles. The difficulty of treating the dynamics of this 
process analytically is well known, and the extensive literature tends to focus on certain 
special regimes where simplifying approximations are possible; for example, instabilities due 
to a cold low density beam ||l| or a gentle "bump on tail" 0-0. 

Early work on the interaction of a narrow spectrum of weakly unstable waves with a 
cold electron beam identified the importance of particle trapping and predicted a relation 
cuf, ~ 7 between the linear growth rate 7 and the bounce frequency u"^ = ekEk/m in the 
nonlinear state that emerges after the linear instability has saturated. [|l| This relation 
is more transparently stated as a property of the electric field E^ of the saturated wave: 
i?fc ~ 7^ as 7 — > 0; a property we refer to as trapping scaling for E}^. Initial studies of the 
saturated state for the bump on tail instability also predicted trapping scaling. [Q,^ These 
investigations all invoke approximations treating the response of the non-resonant electrons 
as linear or adiabatic. 

By contrast, Simon and Rosenbluth constructed a time-asymptotic state by perturba- 
tively expanding the Vlasov equation in the beam density and demanding that secular terms 
vanish at each order. Their procedure led to expressions involving singular functions which 
were defined by prescribing certain regularization procedures; the resulting theory predicted 
nonlinear states with much larger electric fields Ek ~ 1/7. @] Subsequent perturbative cal- 
culations by other groups involve similar assumptions and reach the same basic conclusion 
I^J^ with one exception: Larsen has studied the bump on tail problem using a multiple scale 
expansion in time and velocity; in his formulation, a "singular layer" at the phase velocity 
leads him to posit trapping scaling as an ansatz. |0] Numerical simulations of the instability 
find the initial, possibly metastable, nonhnear state with trapping scaling |P,|| but have 



sometimes claimed to detect a slow growth in Ek on very long time scales. [^] Recent labora- 
tory experiments of an electron beam interacting with an electrostatic wave (supported by a 
travelling wavetube) find the nonlinearly saturated wave amplitude is described by trapping 



scaling over the length scale of the experiment. |T^ 

In this paper we describe a new approach which also treats the Vlasov equation pertur- 
batively, but simplifies the problem in a new way: the initial conditions are restricted so 
that the evolution occurs on the unstable manifold of the equilibrium Fq. 111| Physically 



this restriction means that only the unstable mode is initially excited, rather than viewing 
the mode as one component of an arbitrary fiuctuation. Mathematically the unstable man- 
ifold provides a finite-dimensional setting which partially compensates for the absence of a 
finite-dimensional center manifold in this problem. In the simplest case, the unstable mani- 
fold defines a two-dimensional problem where the invariance of Fq under spatial translation 
implies that the evolution of the mode amplitude decouples from the phase and is described 
by a one-dimensional dynamical system. We formulate and analyze this one-dimensional 
problem. 

The dimensionless equations for the potential 0(x, t) and electron distribution function 
F{x,v,t) are 



oo 



dtF + vd^F + d^(l)d^F = d'j = dvF-l (1) 

J —oo 

where 

/L/2 /-oo 

dx / dvF{x,v,t) = 1. (2) 

-L/2 J-oo 

Given an equilibrium Fq{v, fj,), Eq.(|l|) gives the evolution equation for f{x, v, t) = F{x, v, t) — 

dj = Cf + NU) (3) 

where A/'(/) = —dx<t>dvf and Cf = —vdxf — dx(pdyFQ; here /i denotes parameters such as 
density or temperature that determine the properties of Fq. We assume periodic boundary 
conditions on / thus the Fourier components fk{v) e^^^"^ are discrete multiples of k^ = 27t/L. 



The length of the system is chosen so that as fi varies an instabihty occurs for kc (or 
k = 1) corresponding to an eigenfunction £\I/ = A\& where \E'(a;, v) = e^^"^ ipdv)-, A = —ikcZo, 
and ipc{v) = dvFo/k'^^v — zq). The unstable mode is determined by a root zq = Vp + i'y/kc of 
the dielectric function 

and both Vp, the phase velocity of the mode, and 7 depend on the equilibrium parameters 
/i. For our purposes it is not important exactly what parameters /i represents or how they 
are varied, the mode is assumed to be neutrally stable for fi = fic and to become unstable if 
/i is appropriately shifted away from /i^; the inverse limit ^ —^ ^c from the unstable regime 
will be denoted by 7 — i> 0^. In addition we assume the derivatives d'^ek{zQ)/dz"' have finite 
limits as 7 — > 0+ and that the first derivative e'f^{zo) is non-zero. This latter requirement 
that Zo is a simple root will typically be satisfied when a single parameter is varied, e.g. 
beam density or beam velocity. 

When Fo(f,/i) lacks reflection symmetry as in a beam-plasma system then Vp ^ and 
the eigenvalue A is complex. For 7 > there is a complex conjugate pair of eigenvalues 
(A, A*) in the right half plane and a symmetrically placed pair (—A, —A*) in the left half 
plane; as 7 ^ 0+ this eigenvalue quadruplet merges into the continuous spectrum on the 
imaginary axis. |jl2| For reflection-symmetric problems, Fo{v,fi) = Fo(— f,yu), the eigenvalue 
can be real or complex; for example, two-stream instabilities correspond to real eigenvalues 
of multiplicity two: \E' and \1/* are linearly independent eigenvectors for A. [|13| An instability 
with complex A leads to four-dimensional unstable manifolds when Fq is reflection-symmetric 
and will not be considered here. ||14| 

For A complex (no reflection symmetry) or A real (with reflection symmetry), we decom- 
pose / to isolate the two unstable modes 

/(x, f, t) = [A(t)^(x, v) + cc] + S{x, V, t) (5) 

where A{t) = (\1', /) and (\1/, S") = 0; \1/ is the adjoint eigenvector corresponding to A* and 



{Gi,G2) = J dx J dvGi{x,v)*G2{x,v) denotes the inner product. The equations for A and 
S follow from (^ 

A = XA + i^,J\fif)) (6) 

d^S = CS + Afif) - [i^,Afif))^ + cc] ; (7) 

the linear terms are now decoupled, but nonlinear couplings between A and dfS remain. 
For 7 > 0, the modes \1/ and \I^* span the two-dimensional unstable subspace E^ which is 
invariant under the linear flow dtf = Cf. The nonlinear terms couple the unstable modes 
to the two stable modes (spanning the stable subspace E'^) and to the continuum (which 
spans the infinite-dimensional center subspace E'^); these interactions bend the unstable 
subspace into a two-dimensional unstable manifold which is invariant for the full nonlinear 
evolution. Solutions on this manifold asymptotically approach Fq at an exponential rate e"'^ 
as t -^ — cxD. 

At the equilibrium, the manifold is tangent to E^ and thus can be represented as the 
graph of a function H : E^ ^ E'^ (B E^; see Figure 1. With respect to the decomposition in 
(I), when / is a point on the unstable manifold (denoted f^) then 

r(x, v) = [A-^{x, v) + cc] + H{x, V, A, A*); (8) 

hence the time dependence of S for a solution on the manifold is determined by the dynamics 
of A(t) and the geometry of the manifold: S{x,v,t) = H{x,v, A{t), A*(t)). For such a 
solution dtS = ABaH + A* dA*H and consistency with (j^) requires 

'AdAH + A^dA-H 

Solving this equation for if, then determines the representation of /" in (j^). Setting f = f^ 
in (|^) yields the desired description of the dynamics on the unstable manifold: 

A = XA + i^,Afin) = ViA,A*,fi); (10) 

this is an autonomous two-dimensional dynamical system for A{t). For small A, the graph 
function H is second order in A and /" ~ [A{t)'${x,v) + cc], in this sense the dynamics on 
the unstable manifold corresponds to an initial excitation of only the unstable modes. 



= CH + Af{r)- {^,Af{r))<iJ + cc . (9) 



The spatial translation symmetry of the problem constrains V to have the form 
V{A, A*, yu) = ylp(|y4p, /i) where p is an undetermined function. If Fq is reflection-symmetric 
as in the two-stream instability, then p must be real-valued; otherwise p is complex-valued. 
Finally in amplitude/phase notation A = pe~*^ (pOl) becomes 



p = p Re [p(p^ /i)] p^ = -Im[p(p^p)]; (11) 

thus on the unstable manifold the evolution of the mode amplitude p(t) decouples from the 
phase 6{t) and is a one-dimensional problem. 

For small wave amplitudes, we have investigated the properties of (|ll]) by representing 
p(p^,p) as a power series, 

oo 

p(/,/^) = Ep.(/^)p'^- (12) 

j=0 

Clearly po{fi) = A from (^ and the higher order coefficients Pj(p) are determined by solving 
(y) for if as a power series in {A, A*) then substituting into (0); this calculation will be 
presented elsewhere. ||T5[ Of greatest interest are the properties of pj in the regime of weak 
instability 7^0+: for j > 0, these coefficients are singular at every order 

pAf^) = ^-^[^+om, (13) 

and the remaining asymptotic dependence on the equilibrium Fq is remarkably simple. Each 
bji^Hc) depends on Fq only through the phase e*^'^'^'^^ = efc(fp)*/e^(fp), which is determined 
by the limiting value of efc(-Zo)- More precisely, at each order there is a calculable function 
Qj, independent of Fq, such that 

6,(pJ = g,(e^«('^^)); (14) 

in particular, we find bi{fic) = —1/4 although in general Qj will not be constant at higher 
order. 

The significance of the divergence in (^) is clearer when we introduce a rescaled mode 
amplitude p(t) = 7^r(7t) which varies on the slow time scale t = 'jt and rewrite the 
dynamics ([Tl|) as 
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— = rh+J2Re [b.ifi,) + 0(7)] r^n = ^ - 7 ^ Im [6,(/x,) + 0(7)] r^^' (15) 

where tu = k^Vp is the mode frequency. The 7-^0^ hmit for (|15]) is nonsingular and the 
asymptotic equations have some notable features. First, the dependence on Fq is entirely 
contained in the three parameters 7, uj, and e*^; in particular at 7 = the rescaled amplitude 
dynamics depends only on e*^. If this phase is fixed, then any variations in densities or 
temperatures characterizing Fq do not affect the evolution of r(r). For example, a beam- 
plasma instability (complex A) and a two-stream instability (real A), compared at a fixed 
value of e*^, have identical amplitude equations up to 0(7) corrections. A second feature 
is that the linear term defines a growth rate that is not small, i.e. not 0(7), rather the 
growth rate is unity. Moreover, unless bj happens to vanish, all the higher order terms 
r^-' in the amplitude equation are order unity in the limit 7-^0+. Thus there is no 
small parameter in ([T5|) to justify a truncation of the series, and consequently it is not 
straightforward to calculate the time-asymptotic amplitude r{T ^ 00). However, assuming 
that r{T) approaches a limiting value Too as r ^ 00, then it follows that this time-asymptotic 
state is a BGK mode. [^ 

The implications of (|T^ for the behavior of the electric field \Ei\ at k^ (or k = 1) follow 
from ikcEi = J dv fi{v, t) = A{t) + J dv Hi{v, A, A*); one can show that 

k, \E,it)\ =^\[l+r^ fi(r2, e'«) + 0(7)] (16) 

where Ti represents a power series in r^ whose coefficients depend on Fq only through e*^. 



151 This expression predicts that the trapping scaling \Ei\ ~ 7^ is a universal characteristic 
of the entire evolution; in particular it should hold for the time-asymptotic state. This 
conclusion agrees with the early work on the instability due to a small cold beam, ||I[ and 
generalizes it to a much wider class of instabilities. 

An interesting aspect of (p!5| ) - ([16]) is the apparent absence of the familiar trapping 
oscillations in the evolution of \Ei\. This evolution is determined by r{T) whose dynamics 
is described by an autonomous one-dimensional flow, and it is well known that smooth one- 



dimensional equations cannot describe oscillations. We have conjectured elsewhere [0, on 



the basis of simpler exactly solvable models, that the unstable manifold develops a spiral 
structure away from the equilibrium as illustrated in Figure 2. If this is correct, then 
representing the dynamics on the manifold via a mapping H from the unstable subspace 
yields a vector field on E^ with branch point singularities at the points where the flow 
moves from one branch of the spiral to the next. Thus Re [p(p^,/i)] would have a branch 



point Re [p(p^, /i)] ~ J pi — p^ at pi and as the mode grows p^ir) -^ pi would signal the onset 
of trapping oscillations with the passage of the trajectory to the next branch of the unstable 
manifold. Note that a trajectory will reach such a node in finite time, unlike the more familiar 
situation of a node where the vector field is differentiable and the approach time is infinite. 
In addition the loss of smoothness at p = pb introduces the lack of uniqueness needed by the 
solution to pass through the branch point. Such a spiral structure would present a significant 



obstacle to using the power series (|T2D to determine the time-asymptotic amplitude r^o- 

The coefficients pj (p) are calculated as integrals over velocity and the singular behavior in 
(p^) arises from pinching singularities that develop at f = fp as 7 ^ 0^. The regularization 
procedures proposed in previous treatments would modify these integrals to remove the 
pinching singularity and eliminate the divergences reflected in (JIB]). Any such regularization 
replaces trapping scaling by \Ei\ ~ ^/7 which is typical of a Hopf or pitchfork bifurcation in 
which there is no continuous spectrum on the imaginary axis. In these latter bifurcations the 
coefficients pj are nonsingular, the series for p can be truncated, and time-asymptotic state 
found by balancing the linear term against a cubic nonlinearity. However, our calculations 
show that the singularities in pj found here simply imply a different dependence of p{t) and 
I -El I on 7 as 7 —i> 0^, and that there is no need to introduce an ad hoc regularization. 



Instabilities in other Hamiltonian systems, including ideal shear flows [|T^ and solitary 
waves |T^ , also exhibit key features of this problem, most notably that the unstable modes 
correspond to eigenvalues emerging from a neutral continuum at onset. It will be interesting 
to determine if similar singularities arise in the amplitude equations for the unstable modes in 
these problems. [|l^] By contrast, there is at least one example, a phase model for the onset of 

8 



synchronized behavior in a population of oscillators, in which the critical eigenvalues emerge 
from the continuum at the onset of instability but the amplitude equations are nonsingular 
and ^7 scaling is found (at least in the best understood case of a real eigenvalue). p0|j21 



This difference in the nonlinear behavior seems noteworthy since the linear dynamics of the 
model is qualitatively similar to Vlasov although apparently lacking a Hamiltonian structure. 
22 
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FIGURES 
FIG. 1. Local geometry of the unstable manifold; the equilibrium Fq is at the origin 

FIG. 2. Conjectured spiral structure in the global unstable manifold 
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